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Abstract 

The M a energy which is usually minimized in branched transport 
problems among singular 1-dimensional rectifiable vector measures with 
prescribed divergence is approximated (and convergence is proved) by 
means of a sequence of elliptic energies, defined on more regular vector 
fields. The procedure recalls the Modica-Mortola one for approximating 
the perimeter, and the double-well potential is replaced by a concave 
power. 

1 Introduction 

The name "branched transport" is now often used for addressing all the trans- 
port problems where the cost for a mass m moving on a distance I is proportional 
to I but not to m but sub-additive and typically proportional to a power m a 
(0 < a < 1). The distributions of sources and destinations are given and one 
looks for the path followed by each particle, sums up the mass which moves 
together on each part of the path, and associates to every configuration its 
total cost X)i^ m f- The adjective "branched" in the name stands for one of 
the main features of the optimal solutions: they gather mass together, masses 
tend to move jointly as long as possible, and then they branch towards different 
destinations, thus giving rise to a tree-shaped structure. 

In the case of finite graphs this kind of problems dates back to the '60 in the 
community of operational research (the first paper on the subject is |16|). More 
recently, several different approaches and results concerning the generalization 
to continuous frameworks have been introduced by the community of optimal 
transport. The first paper in this direction is the one by Q. Xia ([21]), which 
will be our main reference. Other approaches have been proposed by Bernot, 
Caselles, Maddalena, Morel and Solimini (see [17l [6] ) . The equivalences between 
the different models as well as a survey on the whole theory are presented in a 
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recent book by Bernot, Caselles and Morel, [7J, who are responsible for most of 
the results therein. 

A satisfactory numerical treatment of these issues is far from being obtained 
and this is one of the reasons for the need for approximation results. This paper 
presents a precise approximation result, in terms of T— convergence (see [H]), 
which can be the origin of possible numerical methods based on PDEs, after an 
idea by Jean-Michel Morel. Actually, the continuous formulation of branched 
transport problems passes through a divergence-constrained formulation (see 
next session, where the model by Xia is quickly sketched) and it is natural to 
approximate it by means of problems which concern more regular vector fields 
(i.e. which are not measures concentrated on one-dimensional graphs, but have 
a density that is at least weakly differentiablc). 

We propose an approximation based on functionals with a concave term and 
a Dirichlet term, with coefficients which increase the weight of the first and 
vanish on the latter as far the approximation parameter e goes to 0. In this way 
one would need to solve a problem of the kind 



under constraints or penalization on V • u. The optimality conditions for this 
problem read as an elliptic vector PDE (a system) with an unknown pressure 
coming from the divergence constraint. In dimension two, the decomposition of 
any vector field into a gradient plus an orthogonal gradient allows to translate 
this into a system of scalar fourth-order elliptic equations. Local minima may 
be found by a gradient descent method passing through a parabolic evolution 
equations, while global ones are more difficult to detect due to the concavity of 
F. Numerical methods for this kind of equations are much more developed and 
efficient rather than the combinatory-based ones which characterize the study 
of finite networks. A partial use of genetics algorithm could thereafter help in 
avoiding local minima. 

Besides possibles numerical applications, the interest of this T— convergence 
result also comes from its comparison with the elliptic approximation of the 
perimeter functional, proposed by Modica and Mortola ([18]) at the beginning 
of T— convergence time. In their case u was scalar and F was a double-well 
potential, enforcing at the limit u (after a suitable rescaling of the functional, 
so that the .F-part has a coefficient going to infinity and the Dirichlet part 
a negligible one) to take values in {0, 1}. And the energy was, in the limit as 
e — > 0, concentrated on a lower-dimensional structure, i.e. the interface between 
the two phases u = and u = 1. The same will happen here: the concave power 
F will play the role of a double well at u — and |u| = oo and the energy will 
concentrate on a one-dimensional graph. 

Being u a vector, in our problem, one could also evoke Ginzburg-Landau 
theory with its approximation (see [H [9] where the problem of the convergence 
for the minimizers of 
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is first addressed). Yet, it is easy to notice, due to the divergence which is a 
bounded measure, that the problem is essentially scalar, as at the limit there 
will be locally one direction only which will be relevant. Actually a singular 
vector measure, concentrated on a lower-dimensional object, must be oriented 
along a tangent direction if we want its divergence to be a measure (otherwise 
it is a first-order distribution). 

Moreover, in dimension two it is possible to take advantage of the usual 
decomposition of a vector field into a gradient plus a rotated gradient so that 
information on the divergence "fix" the gradient part and all the functional 
may be expressed through the other gradient. In this way one would arrive to 
consider the limits of something like 



which are functionals of the form of those studied by Aviles and Giga (Modica- 
Mortola results for higher order energies, see [H [5] where the second order term 
only contains the Laplacian, and lately pQ with the whole Hessian). 

Probably the main goal of the paper is creating a a bridge between two 
different topics in Calculus of Variations: the approximation of free discontinu- 
ity problems on the one hand and the optimization of transport networks. The 
first one, much linked with elliptic PDEs has already been object of high-quality 
researches for decades and is studied in relation with its applications in mate- 
rial sciences and image segmentation. The second is, in its continuous version, 
more recent and linked to the theory of optimal transport by Monge and Kan- 
torovitch, with applications ranging from economics to biology and geophysics 
(see [21]). Elliptic PDEs, dimensional reduction and geometric measure theory 
are very much involved as well. See O [Til HB] for the whole theory and the 
main examples of T— convergence applied to free discontinuity problems: notice 
that here, the discontinuities (or "jumps") are replaced by a bilateral singular- 
ity, in the sense that the rectifiable graphs in the limit problem vanish almost 
everywhere and are concentrated on one-dimensional sets, thus having a double 
jump, whose "intensity" (measured with respect to Tt 1 instead of C d ) enters the 
limit functional. 

The paper will start with two short sections on preliminaries, Section 2 on 
Xia's formulation of branched transport problems, Section 3 on T— convergence. 
Then, we will discuss briefly in Section 4 why to choose the form of the approx- 
imation we will choose, mainly concentrating on the choice of the exponents, 
since they are not obvious (in particular, the function F in (fTJ will be of the 
form .F(|u|) = \u\@ with (3 ^ a). Section 5 will present the detailed proof of 
the main result (r— convergence of the energies, in dimension two only, under 
no divergence conditions). Section 6 will suggest how to apply the result for 
producing interesting approximated variational problems, underlining what can 
be done with the tools we have so far and what could be worthwhile to prove 
in possible future investigations. 
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2 Branched transport via divergence-constrained 
optimization 

We present here the framework of the optimization problem proposed by Xia 
in pMl |2"5] and then studied by many authors (see for instance [7] for a whole 
presentation of the theory). 

Let 51 C M. d be an open set with compact closure SI and A4(fl) the set of 
finite vector measures on S7 with values in M. d and such that their divergence is 
a finite scalar measure, i.e. such that 

sup jy V<f>-du : (f>e C l (U), \\<j>\\ L ~ < lj < +oo (2) 

(as you can see, we do not ask for (f> vanishing at the boundary, i.e. we take into 
account possible parts of V • u on dfl as well) . The value of the supremum in 
f2j) will be denoted by |V • m| (SI), i.e. the mass of the measure "total variation 
of the divergence" of u. On this space we consider the convergence u e — ► u 
u e — u and V • u e — 1 V • u as measures. When a function is considered as an 
element of this space, or a functional space as a subset of it, we always think 
of absolutely continuous measures (with respect to the Lebesgue measure on SI) 
and the functions represent their densities. 

When we take u € M(fl) and we write u = U(M,9,£) we mean that u 
is a rectifiable vector measure (it is the translation in the language of vector 
measures of the concept of rectifiable currents) u = 0£ ■ Ti^ M whose density with 
respect to the Ti 1 — Hausdorff measure on M is given by the real multiplicity 
9 : M — > IR + times the orientation £ : M — > M. d , £ being a measurable vector 
field of unit vectors belonging to the (approximate) tangent space to M at 
Ti. 1 — almost any point. 

For < a < 1, we consider the energy 

M» = {^ Wl ^ = W<U)' (3) 
I +oo otherwise. 

The problem of branched transport amounts to minimizing M a under a 
divergence constraint: 

min{M Q (w) : V • u = /:=/+-/"} . (4) 

The divergence constraint is given in weak form and means 

• du = J 4> d{f--f+) for all G C*°(S1), 

which actually corresponds to Neumann boundary conditions 

V • u = f in SI and u ■ n = on 9S7. 

From now on, we will always think of Neumann boundary conditions when 
speaking about divergences, so that V • u is the linear functional associating to 
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every <f> e C 1 (f2) (independently of the values on 90) the value J V</> • du: if u 
is a regular function this corresponds to a measure which absolutely continuous 
inside f2 with density given by the true divergence, and which has a TC 1 — part 
on the boundary with density given by u ■ n. 

Remark 1. Notice that this is not the original definition by Xia of the Energy 
M a : Xia proposed it in [21] as a relaxation from the case of finite graphs, but 
formula ^ can be seen as a representation formula for the relaxed energy 

M a (u) — inf jliminf E a (G n ) : G n finite graph, u Gn — > u\ , 

where 

E a {G):=Y J <n\e h ), (5) 

h 

for a weighted oriented graph G = (eh, eh, Wh)h (where eh are the edges, e-h 
their orientations, Wh the weights), and uq is the associated vector measure 
given by 

u G := y ^w h e h H\ eh , 

h 

(and the convergence is in the sense of A4(£l)). For the proof of the equivalences 
between the two definition, look at [25] or at Chapter 9 in [7]. 

Notice that in general Problem ([J} admits a solution with finite energy for 
any pair of probability measures (/ + ,/~) (or, more generally, for any pair of 
equal mass finite positive measures), provided a > 1 — 1/d (this is proven in 
|24j by means of an explicit construction). 



3 Variational approximation, preliminaries 

The main result of the paper will be a T— convergence result for a sequence 
of energies approximating M a . We will see in Section 6 that for a complete 
approximation of the problem, one would need to insert the "boundary con- 
ditions" given by the divergence constraints and to prove compactness for a 
suitable sequence of minimizers u e of the approximating problems. 

For precising what we mean by "approximating the energy" and how to use 
the result, let us sketch briefly the main outlines of T— convergence's theory, as 
introduced by De Giorgi (see [13] and [T^]1. 

Definition 3.1. On a metric space X let F n : X — > K U {+00} be a sequence of 
functions. We define the the two lower- semicontinuous functions F~ and F + 
(called r — lim inf and T — limsup F + of this sequence, respectively) by 

F~(x) := infjliminf F n (x n ) 

n — >oo 

F + (x) :— inf {lim sup i^Xn) 

n — >oo 

Should F~ and F + coincide, then we say that F n actually T— converges to 
the common value F = F~ = F + . 
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This means that, when one wants to prove T— convergence of F n towards 
a given functional F, one has actually to prove two distinct facts: first we 
need F~ > F (i.e. we need to prove liminf n F n (x n ) > F(x) for any approxi- 
mating sequence x n — * x; not only, it is sufficient to prove it when F n (x n ) is 
bounded) and then F + < F (i.e. we need to find a sequence x n — > x such that 
limsup„ F n (x n ) < F(x)). 

The definition of T— convergence for a continuous parameter s — > obviously 
passes through the convergence to the same limit for any subsequence e n — ► 0. 

Among the properties of T— convergence we have the following: 

• if there exists a compact set K c X such that mfx F n = w£k F n for any 
n, then F attains its minimum and inf F n — > mini 71 ; 

• if (x n ) n is a sequence of minimizers for F n admitting a subsequence con- 
verging to x, then x minimizes F 

• if F n is a sequence T— converging to F, then F n + G will T— converge to 
F + G for any continuous function G:I->1U {+00}. 



4 Elliptic approximation, intuition and heuris- 
tics 

As we partially mentioned, the result we will present in Section 5 is somehow 
inspired by, or at least recalls most of the results in the elliptic approximation 
of free discontinuity problems (Modica-Mortola, Ginzburg-Landau or Aviles- 
Giga). We will only mention the following (see [18] and [11]) because of its 
simplicity, even if it is probably not the closest one in this two-dimensional 
setting where Aviles-Giga seems closer. 

Theorem 4.1. Define the functional F £ on L 1 ^) through 
F e (u) = 



i / W{u{x))dx + e J\Vu(x)\ 2 dx ifueH 1 ^); 
+00 otherwise. 



Then, if W(Q) = W(l) = and W(t) > for any t ^ 0, 1, the Junctionals F e 
T— converge towards the functional F given by 

J cPer(S) if u = 1 on 5, u = on S c and S is a finite-perimeter set; 
otherwise, 



F{u) 



where the constant c is given by c = 2 y/W(t)dt. 

We precised the value of the constant so that the reader will notice that 
similar constants are involved in our case as well. For the same reason (the 
analogy with the present paper) we precise also that one of the key-ingredient 
in the proof of the above Theorem is the inequality 

1 



-W(u(x)) + e|Vu(a;)| 2 > 2JW {u{x)\\Vu{x)\ = 2|V(i? o u)\, 
e 



G 



where H is the primitive of VW (so that one has F £ (u) > 2TV(H o u)). 
In our study we will instead consider functionals of the form 

E E {u) = e~ yi [ \u(x)fdx + e~< 2 [ \Vu(x)\ 2 dx, (6) 
Jn Jn 

defined on u G K 2 ) and set to +00 outside H 1 C A4(Cl), for two exponents 

71 < < 72. 

As one can see the functional recalls Modica-Mortola's functional to recover 
the perimeter as a limit, where the double- well potential is replaced with a 
concave power. Notice that concave powers, in their minimization, if the average 
value for u is fixed in a region (which is in some sense the meaning of weak 
convergence, i.e. the convergence we use on A4(fl)) 7 prefer either u — or \u\ 
being as large as possible, i.e. there is sort of a double well on zero and infinity. 

We give here a heuristics for determining the exponents (3, 71 and 72. 

Suppose you want to approximate a measure u concentrated on a segment S, 
with multiplicity m, and directed towards one of the direction of the segment, 
via a measure ua whose density is smooth and concentrated on a strip of width 
A around S (for instance by convolution). 

The values of ua will hence be of the order of m/A d ~ 1 and the values of its 
gradient of the order of m/A d . This gives a functional of the order of 

^^-'(^A^'S) 2 - 

In the minimization one will choose the optimal width A, depending on m and 
e, and this gives 

A ~ £2d-/3(d-l) jtj 2d-fi(d-l) • E e ~£ < x > Id- fi(d-l) jyi V z P) 2d-/3(d- 1 ) _ 

The correct choice for a possible convergence result towards the energy ^ which 
is proportional to m a is obtained by imposing that the exponent of m is a and 
the exponent of e is zero, i.e. 

2-2d + 2ad 7l _ (d-l)(a-l) 
~ 3 - d + a(d~ 1)' ~ 2 ~ 3-d + a(d- 1)' 

Notice that 71 and 72 may not both be determined since one can always replace e 
with a power of e, thus changing the single exponents but not their ratio. Notice 
also that the exponent j3 is positive and less than 1 as soon as a G]l — 1/d, 1[, 
which is the usual condition. 

Finally, it is worthwhile to remark that for this choice of exponents the 
dependence of A with respect to e is of the form A ps gfWCd+i) ( an d 72 > 0). 
This implies lim e ^ A = 0, which gives weak convergence of the approximation 
we chose (enlarging u a strip of width A without changing its mass) to u. It 
proves nothing, but it is coherent with the convergence result we want to prove. 

We conclude by underlining the case of dimension 2, since it will be the only 
one we will solve: in this case one has (3 = 4 ^~ 2 and 71/72 = {ex. — !)/(« + !)■ 
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Remark 2. From this heuristics and from the proof that we will present in the 
following, the reader may see that the construction only works for a > 1 — 1/D. 
This is quite astonishing if one thinks that, also for a < 1 — 1/d, there are 
measures which admit possible solution with finite energy (in particular 
when both the measures are supported on a same lower-dimensional set). Yet, 
the problem lies in the kind of approximation we require, which uses measures u £ 
which are more regular and in particular are absolutely continuous with respect 
to C d and may not be concentrated on lower-dimensional sets. 

5 Our main T— convergence result 

This Section will be devoted to the proof of the following theorem: 

Theorem 5.1. Suppose d = 2 and a €]l/2, 1[: then we have T— convergence 
of the Junctionals M™ to cM a , with respect to the convergence of M{$1), as 
e — > 0, where c is a finite and positive constant (the value of c is actually 
c = or 1 (4c a/(l — a)) 1 ~ a , being c = J \ftP — tdt). 

As usual in several T— convergence proofs, we will work separately on the 
two inequalities we need. 

5.1 T — lim inf inequality 

In this part we will consider an arbitrary sequence u e — > u and we will suppose 
that M"(u e ) is bounded. This implies at first that all the u e are H 1 functions, 
and that 

J \u £ f < Ce x ~ a and J \Vu £ \ 2 < Ce" 1 ^ . (7) 
Step 1 Lower bounds on M £ 

Consider a rectangle R C f2. Suppose for simplicity that it is oriented 
according to the x and y axes, i.e. R = [a,b] x [c,d], and set v £ := [(u £ ) x ] + 
(the positive part of the x— component of u £ ) and v' e :— dv £ /dy. Consider 
also a function <f>\ : R — > [0,1], depending on y only, such that <pi = 1 on 
[a, b] x [c + 5, d— S\ and = for y £ {c, d}. For every x, set R x := {x} x [c, d] 
and R' x := {x} x [c + 6, d - 5]. 

Fix a value of a; € [a,b] and let A £ be the maximal value of v £ on R' x 
(which is well-defined for a.e. x) and L £ = A^ 1 , f £ (t) = (t 13 — L £ t) + and 
F £ (t) = J Q f £ (s)ds. One can write 

= f 2 £ {v £ )+L £ v £ -{ V P-L £ v £ )- > f 2 £ {v £ )+L £ v £ ~L £ v £ I RAK = f 2 (v £ )+L £ v £ I K 

where the second inequality comes from the fact that v £ < A £ implies v^—L £ v £ > 
and A £ = m&xn> x v £ . Considering the other term as well one has 

s a -^+e a+1 (v' £ ) 2 ^s^LcvJn, + e a ~ l f 2 {v £ ) + s a+1 (v' e ) 2 

> e a - x L £ v e l R , x + 2e a f £ (v £ )v' £ = e^L.vJ^ + 2e a (F £ (v £ ))'. 
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By multiplying times </>i and integrating on R x with respect to y, one has 

>e a ' x Lj v e dy + 2e a TV(F £ (v £ )cf >1 )-2e a f F^v^dy. (8) 

Since v £ 4>i vanishes at both boundaries and reaches the value F £ (A £ ) inside 
R' x , the total variation 7V(F £ (u £ )0i) is at least 2F S (A £ ) = 2c Q A\ +l3/2 . This 
value may be computed by a change of variable (t = A £ s): 

F £ (A £ ) = [ ° J{tP-L e t)+dt= f \Ja^s0 - L e A e sA e ds = A^+^co 
Jo Jo 

(we used L e = A^ 1 and t 13 > L e t for t < A s ). 

Notice that the last term in ([8]) will tend to zero as e — * 0, after integration 
with respect to x, thanks to Lemma [5.2l below. since (f>\ is a fixed regular function 
and hence <f>' x is bounded, F s (t) < j% s^^ds = ct 1+f} ' 2 and \v £ \ <\u £ \. 

Lemma 5.2. For any bounded energy sequence u e we have 

lime" / M 1+/3/2 = 0. 
Jr 

Proof. First of all write |u e | < 1 + w £ , where w £ = (\u £ \ — 1) , and then 

\u £ \ 1+ ^ 2 <C + C f wl+^ 2 . 



Notice than that |{|u e | > 1}| < § R \u e \P < Ce x ~ a — > 0, thus w £ vanishes on 
a large part of R. This allows to apply standard Sobolev-Poincare inequalities 
HwelU'- < C||w e ||^fi (in dimension two any exponent r < +oo is admitted). 
Remember w £ < \u £ \ and |Vui e | < |Vw e |. 

Now, for any pair of conjugate exponents p and q, we have 



R 



< C^ 1 - a )/P\\i, || 1 +/ 3 / 2 -/ 3 /P < n M-a)/p\\.. nl+P/2-P/p 
— 06 \\ a e\\ L (i+f3/2-f3/p) q Zi N^ellH 1 

< C r £ (l-a)/p-(l+ l S/2- l S/p)( Q +l)/2 > 

Hence we have 

e a [ \u e \ 1+p / 2 < Ce a + Ce~ f *, 
Jr 

where the exponent j p , from the previous computations, is given by 
I- a (2p + /3(p-2))(a + l) a(2-p) 



7 P = a + 



Ap 2p 



It goes to zero provided 7 P > 0, and it is sufficient to choose p < 2 in order to 
get the result. □ 
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Remark 3. Notice that the proof would have been easier if one supposed that 
u £ ■ n = (since one could have directly applied Sobolev-Poincare to u e ) on 
dfl, which is quite natural. Yet, for the sake of generality, we admitted possible 
divergences concentrated on the boundary, i.e.non vanishing values of the normal 
component. Actually, in this H 1 setting, the divergence of u £ is seen as a 
measure on which belongs to L 2 {Q) + L 2 (dft). 

Hence, we will ignore the last term and, for every x, we look at the quantity 

where m' e {x) :— J R , v e dy. To estimate it from below, we will minimize over 
possible values of A £ . We have 

min e a - l A l3 - l m + Ac e a A 1+fi ' 2 = c 2 m a , (9) 

Ae]0,+oo[ 

the minimum being realized by 



\e2c (2 + f3) 



(10) 



We do not precise here the coefficient c 2 appearing in the minimal value but the 
correct computation is the one in the statement of the theorem. Notice that the 
exponent 2/(4-/3) equals (a + l)/3. For computing the last equality in © we 
need to use the relations between (3 and a. 

This is the first point where we start seeing an expression recalling somehow 
M a . 

Step 2 Comparison with M a 

One can call [i £ the positive measure (e Q_1 |u e | /3 + e a+1 \X7u E \ 2 ) ■ C 2 . Since 
d[i e = M"(u £ ), the measures /Lt e stay bounded in the set of positive Radon 
measures on Q. Hence we can suppose /i £ — 1 fx. 

As we did previously, take a rectangle R and keep the same notations. Keep- 
ing the positive x— part only and the derivative with respect to y only one has 
fi £ > (e a - 1 vP + £ a+1 |?4| 2 ) • C 2 , hence we get 

<t>l(vWe> f [m' e (x)] a dx-r £ , (11) 

7? J a 

r £ being a negligible rest corresponding to the last term in ([5]). 

We would like this estimate to pass to the limit as e ~ > 0, using u £ — > u. 
This is not yet possible: to get it, take again a function <p2 ■ R — * [0, 1], similar 
to (pi, depending on y only, but such that <p2 = 1 on [a, b] x [c + 26, d — 26] and 
4>2 (y) = for y e [c, c + 26] U [d — 26, d]. For every x, set 

m' £}+ (x) := / v £ (p 2 dy, m" (x) := / [{u £ ) x ]^(p 2 dy, m £ (x) := / {u £ ) x 4> 2 dy. 

Notice that estimate ([TTT) implies J R (pi(y)dfi £ > f a [m £ + {x)] a dx (up to the neg- 
ligible term r £ ). Moreover, performing for every x the same proof as before with 
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the negative part instead of the positive one, and using max{m" + (x), m" _ (x)} > 
|m e (x)|, one gets 

MvWe> [ \m £ (x)\ a dx. (12) 

R J a 

The measures m e (x)dx are the projections on [a,b] of 02 ■ (w E ) x . They 
obviously weakly converge to the projection of </> 2 u x . Yet, this weak convergence 
is not sufficient for getting the convergence (nor for lower semicontinuity) of the 
right hand side of (fl~2|) . This is due to the non-convex behavior of the function 
Tn i — * |m| Q . To prove this convergence we need more compactness (and hence a 
stronger convergence of (nx)#(<fa(v, e ) x )). 

What we may prove is that the functions x i— > m £ (x) arc uniformly BV in x 
and this will allow for L 1 and pointwise convergence. We will use the fact that 
u £ — > u in A4(f2) implies a bound on V-« £ . Take a function ip : [a, 6] — > K with 
V'(a) = "0(&) = 0. Consider 

m e (x)tp' (x)dx = / (u e ) x (x,y)(f> 2 {y)'ip'(x)dydx 
JR 

= / u e (x,y) ■ V{(j)2{y)^{x))dydx - I {u e ) y {x 1 y)(j)' 2 {y)tl){x) dydx 
Jr. Jr. 

< |v ■ u £ |(n)||^iu» (a , fe) + H^iUiii^iu-^iiV'iu-ct) < c||viUao (ai6) . 

This proves the BV bound we needed and implies that (Tr x )#((j)u x ) is a 
measure on [a, 6] which "belongs to BV" (i.e. is absolutely continuous and has 
a BV density). Let us call m(x) its density. Moreover, one has m £ (x) — > m(x) 
for almost any x. Passing to the limit in (Tl2|) as e — > one gets, by Fatou's 
Lemma, 

4>2(y)dn > / \m(x)\ a dx. 

R J a 

It is quite straightforward that one can let 6 go to and get rid of the 
functions (pi and (f>2 (which actually depend on S). This gives 

fi{[a,b]x]c,d[) > [ \m(x)\ a dx, (13) 



where m(x) is defined as the density at x of the projection tt x of the measure 
u x restrained at [a, 6]x]c, d[. Notice that this measure has a density as one can 
realize by writing it as a series of measures with densities with finite series of L 1 
norms (i.e. we take a sequence of S n going to zero and see that the corresponding 
m n are all L 1 functions, and the series J^n ll m n+i — ^nlli 1 converges since 

En/IOWn+l - {4>2)n\d\u x \ < +00). 

It is useful to see that this estimate implies that u, simply because of its 
attainability as a limit of u £ with bounded M £ energies, is a rectifiable measure. 
This is proven thanks to the following Lemmas 15.31 15.41 

Lemma 5.3. Suppose u £ — > u in M(£l) and M £ (u £ ) < C . Then u is a one- 
dimensional rectifiable vector measure, i.e. it is of the form u — U(M,d,£). 
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Proof. We will use the rectifiability theorem proved through different techniques 
by Federer and White (see [H[23]) and already used to prove rectifiability of 
measures with finite M a — mass (see also [IS]). This theorem roughly states 
that u is rectifiable if and only if almost all its (d — 1)— dimensional slices par- 
allel to the coordinate axis are countable collections of Dirac masses. Take the 
coordinate axis and disintegrate u x with respect to its projection m(x)dx on 
the variable x, thus getting signed measures v x . For every n, divide f£ into 
2 n horizontal strips of equal width and call fi >n (x) the integral of v x on the 
i-th interval strip (without its boundary). The estimate (fT3)) that we already 
established easily gives 

> J m(x) a Y,kn(x) a dx. 

i 

Up to choosing the levels where to put the boundary of the strips, we can ensure 
that for almost any x no mass is given to the boundaries by |u| (if this is not 
the case, simply translate a little bit the strips, and this will happen again for 
a countable set of choices only). Hence, when we pass from n to n + 1, the 
mass of the previous strip is exactly the sum of the mass of the two new strips. 
Notice that, due to \a + b\ a < \a\ a + \b\ a , the sequence n i— > ^ fi, n (x) a is 
increasing. Call G(x) its limit: we have J m(x) a G(x)dx < +00. This implies 
that G is finite almost everywhere. Thanks to Lemma 15.41 almost every v x is 
purely atomic, and the same may be performed on the direction y. This allows 
to apply the white's criterion and proves that u is rectifiable. □ 

Lemma 5.4. For a < 1 and a measure v on a interval (say [0, set 

G£*\v) - sup I E WiknT ; n E N = 

1 +00 otherwise 

G a {v), and in particular if Ga iy) < +00 then v 

Proof. Let us start from proving Ga Up ^ {y) > G a (v): for any n, build a measure 
v n which is purely atomic, with one atom at each point i1~ n and choose the 
mass of such an atom equal to that of v on 7j j71 . This measures converge weakly 
to v, hence we have 

G a (u) < liminf G a {y n ) < G^ up \v). 

Then, we prove the opposite inequality. We can suppose G a {v) < +00. 
Hence v — J2keN a k3x k is purely atomic and for every pair (i,n) we have, 



i i + 1 
2 n ' 2" 



and 

G a {v) = 

Then we have Ga " P ' (v) = 
is purely atomic. 
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thanks to subadditivity, 



W(kn)\ a < £ 



Summing up, we get 

2"-l 



and the proof is obtained by taking the sup over n. 



□ 



Once wc know about u being rectifiable, one can choose rectangles R shrink- 
ing around a tangent segment to the set M at a point x$ (this works for 
Ti 1 — almost any point £0) and get fx> 6 a ■ TC} M , which implies the thesis. 

This kind of proofs follow a standard scheme, i.e. considering measures \x £ 
whose mass gives the value of the approximating energy, and providing estimates 
on the limit measure /i. This estimates are obtained through local inequality on 
J R d/i, so that shrinking R around a point one gets information on the density 
of fi (here it is the density w.r.t. Ti}). Similar proofs (typically comparing /j, to 
C d instead of Ti. 1 ), are quite used in T— convergence problems in the setting of 
transport and location: see for instance [TU] and [2U] 

5.2 T — limsup inequality 

Step 3 The case of a single segment 

Consider the case u = 6 ■ 7~L} S , being S a segment (for simplicity, S = 
[0, 1] x {0}). For producing a "recovery sequence" one can inspire himself at the 
lower bound proof. We will look for a profile u e with the following properties: 

• the x— component only of u £ must be present and must have the same 
sign, so that u £ = v e ei\ 

• Vi> £ = v' e e2 (i.e. v £ only depends on y); 

• the Cauchy- Schwartz inequality used in (jSJ) must be an equality, i.e. one 



• the total variation must actually be given by twice the maximum (i.e. v e 
must be monotone on the two separate intervals before and after reaching 
the maximum: we will realize it by taking a maximal value at y = and 
symmetric monotone profiles around 0); 

• v £ must vanish at the boundary of a certain rectangle, so that one gan 
avoid using the function <pi\ 

• the maximum A £ must be optimal in @; 



needs 
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• as e — > weak convergence to the measure u — 9 ■ H^ s is needed: we will 
realize it by taking different rescaling of the same profile z (say, v e (y) = 
A e z(A e y)), so that A e — ► oo and J R z(t)dt — 9 will be sufficient. 

Look at the conditions that the profile z must satisfy: we need z(0) = 1 and 
z'(t)t < 0, so that the maximum of v E will be A tl realized at zero, and the 
monotonicity conditions as well will be satisfied. Moreover, on t > (for t < 
just symmetrize), z must satisfy 



Thanks to HUJ), the ratio A^ /2 /(eA 2 e ) equals (1 - /3)0/(2c o (2 + 13)). This is very 
good since there is no more e. 

Notice that this equation has a strong non-uniqueness, since the function 
z i ► y?-z is non-Lipschitz. For instance the constant one is a solution, but 
a solution going from 1 to exists as well (just get it by starting from a different 
starting point and see that it has to reach both and 1 in finite time). This 
means that there are several solutions, and J Q z(t)dt may be any possible value 
larger than J °° zo(t)dt (zq being the only solution with no flat part z = 1 around 
t = 0). 

We want now to check that 9 is larger than the lower bound for the integrals, 
and this gives the final property that z needed to satisfy. 
We need 

9 f°° . . , f 1 z (1-0)6 

- > / z (t)dt = / v — — dz 

2 - Jo V ' Jo V^I(2 + /3)2c 

(the integral has been computed by change of variables z = zo(t), which is 
possible for the solution zo only, since it is injective). This means that we 
must compare Co := J —j==dz and Co = J \fzfi — zdz and prove Co < 

c (2 + /3)/(l-/3). 

Compute Co by integrating by part: 

f\ n — a f 1 ^ - 1 a Q—M r 13 

co = 1 • V z p - zdz = - I z — dz = — - — Go - -co, 

Jo Jo 2v zP — z £ ^ 

which implies Co = Co ^§ an d is enough to get the desired inequality. 

Obviously, one needs after that to perform a correction of u s near x = 
and x — 1 so that the function actually belongs to H 1 (in order to avoid 
discontinuities at the two ends of the segment), and to control the extra energy 
one pays, as well as the divergence. 

One possibility for the case of the single segment is the following: the profile 
we got for v £ is of the form A e z(A e t) and one can simply replicate it radially on 
a half disk, so as to ensure regularity. In this case we need to estimate. Denote 
for simplicity by / E (r) the radial profile one performs and by B £ the half ball 
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where it is non-null. We need to estimate 

e*- 1 I \fef + s a+1 [ \ff, [ \f' E \, and / |/ e | 

JB E J B E JB C JB e 

Just use |/ e | < As, |/ e | < CA\ and \B e \ < CAj 2 (and the value of A e which 
is of the order of e~ 2 ^ 4 ^^): both terms in the energy will be of the order of 
e ( Q +!)/3 ; fa e second term to be estimated (the divergence) will be bounded and 
the last will be of the order of e( Q + 1 '/3 as we u_ 

In this way we have produced a sequence u e that converges to u (weak 
convergence is ensured by construction, weak convergence of V • u e comes from 
the bound we just proved). 

Step 4 Conclusion by density 

As we pointed out in Section 2, it is well known from the works by Q. Xia 
on (see [24]) that the energy M a is obtained as a relaxation of the same energy 
defined on finite graphs. This implies that the class of finite graphs is "dense 
in energy" in the space A4(fl) and T— convergence theory guarantees that it is 
enough to build recovery sequences for such a class (see [H]). One can also 
impose the condition GO dil = 0, so to avoid problems at the boundary. 

For dealing with u = uq, where G is a finite graph, one can simply consider 
separately the segments composing G and apply the previous construction of the 
previous step. Possible superpositions of the part of u £ coming from different 
segments will happen only on regions whose size is of the order of A~ Y and hence 
negligible in the limit (and on such a region, one can use |ui+U2|' 3 < I'l^l^ + I^l' 3 
and |V(ui+U2)| 2 < 2| Vui | 2 + 2| Vui | 2 ). Not only, the number of nodes will be 
finite and hence the bound on the divergence will stay valid. 

Improvement Better connections at the junctions 

The construction for the T — lim sup that we just detailed provides a sequence 
u e converging to u in M. (Q) which works very well in the case of a single segment 
but which could be improved in general. Actually, in the case of a single segment 
the divergence of u e has the same mass as that of u, since it replaces two Dirac 
masses at the two extremal points with two diffuse masses, concentrated on half- 
balls of radius cA~ Y around the points. The mass of u £ itself is only slightly 
larger of the mass of u, due to the part we added at the extremities. 

This changes a lot when one considers more than one segment, since the 
divergence of u £ at the nodes will be given by the superposition of different 
densities, each corresponding to a segment: the integral will vanish due to 
compensations, but ||V • UelU 1 could have increased a lot. 

This may be corrected thanks to the following lemma. 

Lemma 5.5. If g G C 1 {Br) is a function with zero mean on a ball of radius 
R, then there exists a vector field in C 1 (Br) such that V • v — g, v = on 8Br 
and 

IMU- < C {R\\g\\L- + i? 2 ||V 5 ||^) ; ||V«||| a < C (i? 2 |MII~ + R^gWl^) , 
where C is a universal constant. 
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Proof. Take w the solution of the elliptic problem 



I Aw = g in B R , 
[ w = on 3-Br. 

Thanks to J g = one has J aB dw/dn — 0. This allows to define a function 
4> : dBft — ► K with <fi' = du/dn. Take then a cut-off function x(r) such that 
X(r) = 1 if r € [2R/3,R] and x (r) = if r e [0,72/3], | X | < 1, Ix'l < C*/i?, 
|x"| < C/R 2 and define ^(a) = x{\x\)4>{ Rx l\ x \)- 

Now take u = Vw + Rot(Vip), where i?o£ denotes a 90° clockwise rotation. 
In this way V-v = V- Vw = g, since the rotated gradient part is divergence-free, 
and both the normal and the tangential component of v on the boundary vanish 
(since the tangential component of V?/> compensates the normal one of V and 
the tangential component of \7w and the normal of V"0 are zero) . 

We only need to check the bounds on the norms. These bounds come from 
standard elliptic regularity theory (see for instance [T5]), since one has 



\Vw\ 




< 


C(R\\g\\ L -- 


f R 2 \\Vg\\L~) 


1101 




< 


Ci?||V0|| L ~ 


< CR\\Vw\\ L 


D 2 w 


\\h 


< 


C\\Aw\\l 2 = 


C\\g\\h, 


\D^ 


\\h 


< 


CR\\D 2 w\\h 


(OB R ) < CR2 \ 






< 


CR 2 \\Vg\\ 2 L , 


+ C\\g\\lz. 



The last line of inequalities come from the combination of a trace inequality 
in Sobolev spaces applied to D 2 w (where the two coefficients of the gradient 
and the function part have different scaling with respect to R) with a regularity 
estimate for Dirichlet problems (estimating the H k+2 norm of the solution with 
the H k norm of the datum). Combining al the ingredients give the desired 
estimate (we pass to the L°° norms for the sake of simplicity). □ 

With this lemma in mind, we can notice that at every node, the divergence 
of the function u £ we gave before is composed by the zero-mean sum of some 
functions <?i of the form f'(r)(x ■ ei)/r (this is the divergence of the vector field 
directed as e% with radial intensity we used above), where e, is the direction 
of the corresponding segment. Each function is supported on a half ball whose 
simmetry axes follows and the radius is of the order of A~ x w e 2 ^ 4 ^^. 
Setting for each node g = J^i 9i one nas J 9 = 0, the support of g is included in 
B R with R m e 2 /(4-/3) an d \ \ g \\ Lx , < CRT 2 and ||Vff||i,- < CR~ 3 . 

Hence, one can add at every node a vector field v as in Lemma 15.51 so as 
to erase the extra divergence. This would cost no more than something of the 
order of 

e a ' l R 2 \ M |£» + e Q+1 1 |Vw||| 2 < Cs^R 2 -? + Ce a+1 ± = Ce^'\ 

Also the mass of the vector field is not that changed, since we only added 
R 2 \\v\\l°° <C £ ( q+1 >/ 3 . 
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6 Applications and perspectives 



Besides the interesting comparison aspects of this result with respect to the 
similar ones in the approximation of free discontinuity problems, one of the main 
goal of this study concerned possible numerical applications, as we mentioned 
in the introduction. 

We want to replace the problem of minimizing M a under divergence con- 
straints with a simpler problem, i.e. minimizing M". 

The idea would be solving 

min{M«(u) : V • u = f s } , (14) 

being f £ a suitable approximation of / = / + — / — ., and proving that the mini- 
mizers of (fl4|) converge to the minimizers of ((4]) . 

Theorem 15.11 proves a T— convergence result, which should give the con- 
vergence of the minimizers, but the problem is that we never addressed the 
condition V • u = f e , nor we discussed the choice of f e . We will come back to 
this question, which is still open, later. 

What we will consider here are penalization methods, which are quite natural 
as far as numerics is concerned. One could decide to replace the (quite severe) 
condition V-u = / at the limit with a weaker one, concerning a distance between 
V • u and /, possibly with a very high penalization coefficient. Considering 
problems of the form 

min M?(u) + G(V-u;f), (15) 

where G is a functional defined on pairs of finite measures, which is continuous 
w.r.t. weak convergence, and such that G((fi, v), (f + ,f~)) — if and only if fj, = 
v, would play the game, since the T— convergence result could be easily extended, 
thanks to the last property of this kind of convergence that we presented in 
Section [3] 

Yet, this would converge to 

min M a (u) + G(V-u;f), (16) 

which in general is not exactly the same as imposing V • u = f. If on the 
contrary one wants a penalization of the kind £ G, so that at the limit one 
needs G = 0, i.e. V • u = /, then the whole T— convergence result would have 
to be re-established (and would present some difficulties), since we are no more 
in the same framework of adding a same continuous functional. 

Hence, we will stick to the case of a single penalization and we will start from 
the following consideration which works quite well if the divergence is given by 
the difference of two fixed mass positive measures. 

As explained in [24j . the quantity d a (fj,, u) := min{M a (it) : V-it = fi — v\ 
is a distance on the set of positive measures with the same finite mass that 
metrizes weak convergence. 

As a consequence of d a being a distance, due to triangular inequality, solving 

min _ 2d a {f + ,n + ) +min{M Q (u) : V ■ u = /i+ - + 2d a (f~, pT) 
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amount to choosing /u ± = f ± and solving min{M a (u) : V • u = /+ — /"}. 

This means that choosing 2d a as a penalization is a clever strategy (even if 
it only allows for considering probabilities or, in general, hxed mass measures). 

To describe a precise problem where we can insert this penalization we need 
to introduce a slightly different space. Consider the space Y(il) c x 
V(U) x V(Q) defined by 

Y(Q) := {(u, /i, v) e M(Q) x V(Tt) x V(Q) : V • u = [i - u}, 

endowed with the obvious topology of componentwise weak convergence. 
On this space consider the functionals 

(u, fx, u) i ^ M°(u) + 2d a (f + , /J,) + 2d a (/- )M ). (17) 

It is an easy consequence of our previous T— convergence result and of the fact 
that the quantity G((fi, v), (/+, /")) := 2d Q ,(/+, / u) + 2d a (f~,(i) is continu- 
ous on this space that we still have T— convergence to the functional M a (u) + 
G((/Lt, v), (/ + , / _ )). The minimization of this last functional being equivalent 
to the minimization of M a (u) under the divergence constraint V ■ u = fi — v one 
has obtained a useful approximation of any branched transport problem. 

Notice that in this framework the imprvement that we proposed at the end 
of Section 5 is crucial: without this kind of correction on the divergence, one 
would got approximating measures u e (in the T — lim sup construction) whose 
divergence is no more the difference of two probabilities, since some extra mass 
at every node is added both to the positive and to the negative part. 

Nevertheless, the approximation is not satisfactory yet, especially for numer- 
ical purposes: we are trying to suggest methods to approximate the minimiza- 
tion of M a and we propose to use the distance d a itself!! For computing it one 
should probably solve a problem of the same kind and no progress would have 
been done. 

A possible escape we suggest is replacing d a with other quantities, which 
are larger but still vanish when the two measures coincide. This is possible 
for instance thanks to [19], where the inequality d a < CW-j 2 " -1 is proven, W\ 
being the usual Wasserstein distance on V(Pl). We will call G\ the functional 
Gi((n,v), (/+,/-)) := CW^ a - 1 (^f+) +CW? a - 1 (vJ-). The distance W x 
has the advantages of being independent of the branched transport problem 
and of admitting more numerical methods for computing it (it is, by the way, 
a dual Lipschitz norm). Notice that other possible attachment terms, easier to 
compute, such as the L 2 norm, are not adequate because they are not continuous 
with respect to the convergence we have (and because in general we do not want 
to require / G L 2 ). 

It is anyway important to stress that, the approach on without penal- 

ization stays useful for a lot of problems where the divergence is not prescribed 
but enters the optimization (think at min M d a (fi, v) + F(ii)). Some of this prob- 
lems are addressed in |22) . for instance for urban planning or biological shape 
optimization. 

Yet, even if there are, as we showed, was to overcome the problem, we find 
anyway interesting to ask the following question: 
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Open question 1. Given f , is it possible to find a suitable sequence f s f 
so that one can prove T— convergence of the junctionals u i— » M"(u) + 7v-«=/ e 
to u i > M a {u) + Iv-u=f (being I the indicator function in the convex analysis 
sense, i.e. +00 if the condition is not verified or zero if it is)? is it possible to 
find f e explicitly, for instance as a convolution of f with a given kernel? 

The second issue we want to address, after the one concerning divergence 
constraints, deals with the convergence of the minimizers. T— convergence is 
quite useless if we cannot deduce that the minimizers u £ converge, at least up 
to subsequences, to a minimizer u. Yet, this requires a little bit of compactness. 
The compactness we need is compactness in M(Q), i.e. we want bounds on the 
mass of V • u e and of u e . The first bound, has been guaranteed by the fact that 
we decided to stick to the case of difference of probability measures. On the 
contrary, the bound on |ti e |(fi) has to be proven. 

Notice that M"(u £ ) < C will not be sufficient for such a bound, as one 
can guess looking at the limit functional. Think at a finite graph with a circle 
of length I and mass m on it: its energy is m a l which provides no bound on 
ml (its mass), if m is allowed to be large. Actually, what happens on the 
limit functional is "bounded energy configuration have not necessarily bounded 
mass, but optimal configuration do have" . This is due to the fact that, if / + and 
/~ are probabilities, then m < 1 on optimal configurations (and no cycle are 
possible, by the way). Notice that this statement does not depend on m m a 
being concave, but simply increasing in m. The same kind of behavior is likely 
to be true on the approximating problems, but after months of work a proof of 
this fact did not appear. 

Hence we conclude with two points: we ask a question and we give a final 
suggestion for getting a useful T-convergence result. 

The suggestion is quite naive : just take a sufficiently large number K so that 
every minimizer u of the limit problem satisfies H(f2) < K and then restrict 
the analysis to the compact subset Yk(£Y) := {(u,fi,v) € Y(il) : \u\(£Y) < 
K}. Unfortunately, there exists no continuous functional *f> on F(f2) which is 
coercive on the mass of u. If it was the case, one could have added a term f 
the form F(u) = (^(u) — K) 2 + to the limit and approximating functionals. It 
is on the contrary a reasonable question, simpler than the one we are going to 
detail below, whether there exists a continuous functional VE' such that, for the 
minimizers u e , the condition ^(u e ) < C implies a bound on the mass (taking 
advantage of possible extra estimates on u £ ). 

Thus, the problems 

min {M £ » + Gi((/i, u)(f+, /")) («, ft, v) e Y K (Q)} 

approximate (with T— convergence) the limit problem given by miny^Q) M a (u)+ 
Gi((/z, v), (/+, /")), which is equivalent to min M a (u) : V • u = f. 

Here as well a proper use of the estimate given in the last improvement in 
Section 5 are important. Actually, one needs to produce a recovery sequence 
u £ for the T — lim sup construction which does not violate the mass constraint 
given by K. This may be obtained in the following way: instead of considering 
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as a dense subset the set of all vector measures u corresponding to finite graphs, 
take the set of all such measures u satisfying \u\ < K. This set is still dense 
in energy. Afterwards, one notices that the construction which was provided 
satisfied |u £ |(0) < + Ce a+1 (where the constant C depends actually on 

u, but not on e, since it depends on the number of nodes and of branches at 
every node). This means that, if u satisfies the strict inequality \u\ < K, this 
will be true for u e as well for e sufficiently small. 
And, as obvious, the question is 

Open question 2. Prove a bound on the L 1 norm of the minimizers u e (or 
on suitable minimizers u e , if needed). If possible, prove it for minimizers u E 
which minimize under a divergence constraint V • u £ = f s (so that it will 
be true even if we add penalizations on the divergence). Is is possible to do it 
simply by writing and exploiting optimality conditions under the form of elliptic 
(second or fourth order) PDEs? do we need to pass through more geometrical 
tools such has some surrogate of the no-cycle condition? Is this bound equivalent 
to a bound on some m a (being m the flux thrugh some segments or curves, as 
in what we proved in Step 1, where m(x) is the flux through R x ) coupled with 
some bound on the fluxes m? 

For finishing this questioning section, here is the last natural one: 

Open question 3. Prove the same results as in this paper or investigate what 
happens in M. d , for d > 3. 

Notice that our proofs almost never used (up to the last improvement in 
Section 5) any gradient decomposition of vector fields in R 2 such as u = V<f> + 
RotVip. Yet, this would have been a typical trick for managing divergences in 
two dimensions. In our opinion the point where we used the most the fact that 
we are in M 2 is when we disintegrate with respect to x and we estimate the total 
variation in y by the oscillation (which is very one-dimensional). 
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